Ordering in Cs 2 CuCl4 : Is there a proximate spin liquid? 
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The layered spiral magnet Cs2CuCk displays several interesting properties that have been sug- 
gested as evidence of proximity to a two dimensional quantum spin liquid. In this paper we study 
a concrete version of this proposal and suggest experiments that can potentially confirm it. We 
study universal critical properties of two-dimensional frustrated quantum magnets near the quan- 
tum phase transition between a spiral magnetic state and a spin liquid state with gapped bosonic 
spinons in the framework of the 0(4) invariant critical theory proposed earlies-. Direct numerical 
calculation of the anomalous exponent in spin correlations shows that the critical scattering has 
broad continua qualitatively similar to experiment. More remarkably we show that the enlarged 
0(4) symmetry leads to the same slow power law decay for the vector spin chirality and the Neel 
correlations. We show how this may be observed through polarized neutron scattering experiments. 
A number of other less dramatic consequences of the critical theory are outlined as well. 
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I. INTRODUCTION 



Despite the great amount of progress made in our the- 
oretical understanding of spin liquid phases^** of quan- 
tum magnets, there are to date no convincing experimen- 
tal examples of such phases in spatial dimension d > 2. 
A particularly promising candidate - much discussed re- 
cently - is the material CsiCuCl^. This is a layered 
Mott insulator whose spin physics is accurately mod- 
elled by a spin- 1/2 Heisenberg antiferromagnet on an 
anisotropic triangular lattice. Neutron scattering mea- 
surements have revealed several unusual properties of the 
spin excitations* which have motivated intensive theo- 
retical research works§i£*2*i£. There is incommensurate 
spiral magnetic order at low temperatures (below 0.62K). 
The low energy excitations of this ordered state are spin 
waves. However spin wave theory fails to quantitatively 
describe the spin wave line width seen in inelastic neutron 
scattering. In particular there are long tails that extend 
to reasonably high frequencies in the inelastic scattering 
spectrum. This structure may be understood as indica- 
tive of a broad continuum of excited states. The con- 
tinuum survives upon heating above the magnetic phase 
transition (though the spin wave itself does not). In ad- 
dition the phase diagram in an external magnetic field 
shows several interesting features. 

The broad continuum is reminiscent of spinon exci- 
tations in one-dimensional Heisenberg antifcrromagnetic 
chain. However the continuum disperses in both spatial 
directions so that a two dimensional description of the 
physics may be more appropriate. A number of workers 
have therefore suggested interpretations of the neutron 
data in terms of fractional spin (spinon) excitations in 
two dimensions into which the magnons decay. In par- 
ticular Coldea et alA suggested that the material may 
be viewed as being close to a quantum phase transition 
between a spiral Neel state and a spin liquid state. The 
purpose of the present paper is to explore this possibil- 



ity in greater detail and to suggest concrete experimental 
tests to confirm (or rule out) this proposal. 

We begin with some general considerations. It is by 
now established that a number of different kinds of spin 
liquid phases are theoretically possible in two dimen- 
sional quantum magnets. We will restrict attention to 
a particular kind of spin liquid phase whose excitations 
consist of bosonic spin- 1/2 spinons with a full spin gap. 
In addition there are gapped Zi vortices (visons) which 
act as sources of it flux for the spinons. As shown many 
years ago in Ref. Q] such a spin liquid phase admits a di- 
rect second order transition to the spiral Neel state which 
is simply driven by condensation of the bosonic spinons. 
Specifically the transition was argued to be in the uni- 
versality class of the classical (0(4) fixed point in three 
dimensions. Here we point out several remarkable conse- 
quences of this theory which may be used to test its ap- 
plicability to CS2CUCI4 or other materials. Perhaps most 
interestingly we show that the large extra 0(4) symmetry 
that emerges at the critical fixed point unifies seemingly 
different competing orders. The Neel vector correlations 
have the same slow power law decay as the vector spin 
chirality. This non-trivial prediction can potentially be 
probed in polarized neutron scattering experiments. 

In passing we note that other more exotic spin liquid 
phases with (possibly gapless) fermionic spinons could 
potentially exist as stable phasesii*i2iiii^ but the tran- 
sitions to the Neel state have not been studied. The 
corresponding critical theories are also likely to be more 
exotic - we will therefore defer consideration of such spin 
liquid phases and the corresponding critical points for the 
future and focus here on the simpler case with gapped 
bosonic spinons. 

A controlled calculation in which both the magneti- 
cally ordered state and such a gapped spin liquid state 
appear is provided by considering a large- N generaliza- 
tion of the S =1/2 Heisenberg modeli on the anisotropic 
triangular lattice. Specifically, the large-A^ limit of a 
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bosonic Sp(7V) Heisenberg modeli£*i& was used to ob- 
tain the mean field phase diagram as a function of the 
ratio of intra-plane anisotropic couplings, J'/ J, and the 
strength of quantum fluctuations, re = "2S"', where re 
plays the same role as 2S (= value of spin at each site) 
in the SU(2) limit. It was found that, for a range of 
J'/ J, the semiclassical limit (large re) leads to a spiral 
magnetic order. As the quantum fluctuations become 
stronger (small re), a spin liquid state (Z^ spin liquid) 
with gapped deconfined bosonic spinons emerges. 

The possibility that CS2CUCI4 is proximate to the 
quantum critical point between the incommensurate spi- 
ral and spin liquid states will affect the behavior at in- 
termediate energy/length scales. As mentioned above 
Chubukov et aim showed quite generally that the transi- 
tion was in the universality class of the 0(4) fixed point in 
three euclidean dimensions. This is much higher symme- 
try than in the original microscopic model. This enlarged 
0(4) symmetry acts naturally on spinon degrees of free- 
dom which thus emerge as the useful variables already 
at the transition to the spin liquid. The theory predicts 
a large anomalous exponent fj for the spin-spin corre- 
lation function since the spins are composite operators 
in terms of spinon operators. Such composite operators 
usually have large anomalous dimensions. A qualitative 
physical picture is simply that the spin-1 magnons decay 
rapidly into the spinons thereby leading to broad line- 
shapes in the inelastic neutron scattering. However right 
at the critical point the spinons are not free particles. 
From the large- N calculation extrapolated to the phys- 
ical case of SU(2) spins, the anomalous exponent may 
be estimated to be 1.54m. A more accurate value of fj 
can be obtained directly by classical Monte Carlo sim- 
ulations of the 0(4) nonlinear sigma model in three di- 
mensions. In the present work and in Ref. Il7l it is found 
that fj = 1.373. It is worth noting that the neutron 
experiments^ in CS2CUCI4 were fit to functional forms 
that also suggest large anomalous dimension 77E in the 
range 0.7 to 1. There is however reason to question these 
measurements of the actual numerical value of fj as we 
discuss in Sec E] Nevertheless the large anomalous di- 
mension predicted by the theory is qualitatively consis- 
tent with the broad lineshapes seen in experiment. 

Thus it is certainly desirable to have other qualitatively 
distinct predictions of the critical theory. It is our hope 
that the enlarged 0(4) symmetry and the consequent en- 
hanced vector spin chirality correlations will provide such 
a sharp test. We also outline some other less dramatic 
consequences of the extra 0{4) symmetry that too may 
be useful. 

From a theoretical point of view our considerations 
are rather similar to those in a recent studyA& of "alge- 
braic spin liquid" phases where too the low energy the- 
ory is characterized by non-trivial enlarged symmetry as 
compared to the microscopic model. There as in the 
present problem this enlarged symmetry acts naturally 
on "spinon" degrees of freedom and leads to non-trivial 
relationships between the fluctuations of rather different 



competing orders. Perhaps such phenomena are common 
in correlated systems. 

The rest of this paper is organized as follows. In 
Sec. [n] we introduce the model that we use to describe 
CS2CUCI4. We briefly review the 0(4) invariant critical 
theory in Sec. lIIII and discuss the anomalous exponent in 
Sec. HVAl In Sec lTVBl we study the vector spin-chirality 
and scalar-chirality correlation functions. In Sec we 
discuss possible experiments to test our predictions and 
the critical theory. Finally, we summarize our results and 
conclude in Sec IVII 

II. MODEL 

The magnetic ions Cu 2+ in Cs 2 CuCl4 carry 5=1/2 
spin moments that reside on a stack of triangular layers^. 
The intralayer antiferromagnetic interaction (the interac- 
tion in the b-c plane; see Fig.l) is anisotropic with two 
coupling constants: J ~ 0.375meV along the b direction 
and J' ss J/3 along the zig-zag bonds^ (see Fig.^l. The 
interlayer interaction (the interaction along the a direc- 
tion) is weak£ with a coupling constant J" = 0.045 J <C 
,/. Therefore the system is quasi two-dimensional and 
can be modeled by a two-dimensional frustrated Heisen- 
berg Hamiltonian on the triangular lattice. The main 
part of the Hamiltonian reads 

H = J^8i-8j + J' -S.- (1) 

(id) {(id)) 

Here Si are spin- 1/2 operators at the sites i of a two 
dimensional triangular lattice. The first sum runs over 
bonds in the b direction and the second sum runs over 
the zig-zag bonds. There is also a weak Dzyaloshinskii- 
Moriya interaction^. We will neglect this term when we 
consider the critical theory and discuss later its role in 
experiments. 




FIG. 1: Schematic representation of the triangular lattice 
in the b-c plane. The a direction is perpendicular to the 
plane. There are two coupling constants; J along the hor- 
izontal chains and J' along the zig-zag bonds. 

The weak interlayer interaction stabilizes magnetic 
long range order below T = 0.62K 20 . The order is 
incommensurate because of the frustrated anisotropic 
interactions in the triangular planes and occurs at an 
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incommensurate wave vector Q = (0.5 + eo)h* with 
e = 0.030(2) and b* = (2^/6,0,0). The correspond- 
ing classical ordering wavevector is Q = (0.5 + e c )b* 
with e c = (1/tt) arcsin(J'/2J) = 0.053. The substan- 
tial difference is due to large quantum fluctuations that 
renormalize the ordering wave vector. 

A. Order parameter 

The spiral ordering pattern may be represented using 
two orthonormal unit vectors ri! and n 2 through: 

(Sfa)) oc ni cos(Q • n) + n 2 sin(Q • r;), (2) 

where tii and n 2 satisfy 

m • n 2 = 0, n\ = n\ = 1. (3) 

The 111,2 are clearly vectors under global spin rota- 
tion, and together they define the order parameter for 
the spiral state. In addition to global spin rotations the 
spin Hamiltonian Eq. ^ is invariant under various lattice 
symmetries and time reversal. It is useful and straight 
forward to work out the transformation properties of ni j2 
under these symmetries. We find 

Tj : (m + m 2 ) -» (m + m 2 )e ,;Q! , 
R C ,I : (m +in 2 ) (ni - in 2 ), 

T : (m +in 2 ) -(ni - in 2 ). (4) 

Here Tj is a unit translation along a lattice vector 1 direc- 
tion, R c is reflection about the c axis, I is lattice inver- 
sion, and T is time reversal. Other symmetry operations 
may be obtained as combinations of the above. 



III. THEORY OF QUANTUM TRANSITION 

n 

In this section we briefly review the theory of Ref. 1 
for the quantum transition between the Neel and spin 
liquid states. We will also use this as an opportunity to 
clarify several possible confusions with other superficially 
similar results in the literature. 

The set m, n 2 defines the order parameter for the spi- 
ral state. To study phase transitions out of the spiral 
state it is necessary to allow fluctuations where 111,2 vary 
slowly as a function of space and time. The symme- 
try properties in Eq. 21 actually imply that the resulting 
continuum theory has an extra U(l) symmetry over and 
above that of spin rotations. This may be seen by con- 
sidering translations along the b direction. We have 

(ni + in 2 ) - (ni + in 2 )e imQ * (5) 

for a translation by m lattice sites. With incommensurate 
Q x , this effectively translates into a full U(l) symmetry 
that rotates between the ni and n 2 fields. 



The order parameter manifold defined by ni )2 al- 
lows for topological vortex defects with a discrete Z 2 
character 21 . These Z 2 vortices have energy logarithmic 
in the system size in the ordered state due to the long 
distance distortion of the order parameter. Now consider 
disordering the spiral order while keeping the core en- 
ergy of these Z 2 vortices finite. Ref. U argues that the 
resulting state is a fractionalized spin liquid with bosonic 
spin-1/2 spinons. The Z 2 vortices survive into the spin 
liquid phase but now only cost finite (not divergent) en- 
ergy. 

A theory for this transition is obtained by writing 

ni + in 2 = e Q7 z 7 <T Q/3 z /3 , (6) 

where e is the antisymmetric tensor, <j a are the Pauli 
matrices, and z a = (z-\,zi) a two component complex 
unit vector satisfying 

zJz = 1. (7) 

It is easy to check that the parametrization © satis- 
fies the constraints given by Eq. |3| The z a transform as 
spinors under the SU (2) spin rotation and describe spin- 
1/2 spinons. This representation clearly has a Z 2 gauge 
redundancy associated with changing the sign of z at any 
point in space-time. Thus a reformulation in terms of z 
may be fruitfully viewed as a theory of the z fields cou- 
pled to a Z 2 gauge field. The corresponding Z 2 gauge 
flux is associated with the Z 2 vortices discussed above. 
These vortices stay gapped at the transition between the 
Neel and spin liquid states and may hence be ignored 
for a low energy description of the critical point. Thus 
we may obtain the critical behavior by focusing only on 
the spinons and ignoring their coupling to the Z 2 gauge 
field. Detailed arguments show that the critical theory is 
in fact in the 0(4) universality class in D = 2 + 1 dimen- 
sions where there is full rotational symmetry between the 
four real numbers described by (zf , z{). In effect the ex- 
tra (7(1) symmetry of rotation between ni and n 2 has 
been enlarged to SU(2). Combined with the SU(2) spin 
rotations the full symmetry is 577(2) x SU{2) ~ 0(4). 
Thus the critical properties may be computed using the 
Euclidean action 

S= f d 2 xdr ]T -dnzldpz*. (8) 

The relationship between the z and ni, 2 fields may also 
be expressed in a different way that will be fruitful later. 
Let us introduce an SU(2) matrix U built out of z: 

Then U satisfies 

U^a a U = R ab a\ (10) 
where R is a 3 x 3 rotation matrix. Clearly 

R ab = \tr (U^a a Uo b ) . (11) 
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It is readily checked that ni and n 2 are the first and 
second column respectively of this rotation matrix. The 
third column represents the unit vector 

n 3 = ri! x n 2 (12) 

which is orthogonal to both ni and n 2 . Symmetry under 
physical spin rotations correspond to left multiplication 
of U by an SU (2) spin rotation matrix: 

U -> VU (13) 

with V E SU(2). The enlarged 0(4) symmetry implies 
that right multiplication 

U^UV (14) 

is also a symmetry of the critical fixed point. 

Several comments are in order on these results. Prior 
to Ref.Q, analysis of a continuum non- linear sigma model 
appropriate for non-collinear magnets had led to the sug- 
gestion of enlarged 0(4) symmetry in 2 + e space-time 
dimensions 22 . However numerical calculations on stacked 
triangular lattices near their finite temperature ordering 
transition failed to observe the predicted 0(4) universal- 
ity class2£. Rather the evidence sup ports a transition in 
the "chiral" universality class of Ref. 23 24. How then are 
we to reconcile this with our claims about the quantum 
Neel - spin liquid transition being in the 0(4) universality 
class? The answer lies in the nature of the paramagnetic 
phase. In the classical stacked triangular lattice the natu- 
ral paramagnetic phase is the trivial one that is smoothly 
connected to the high temperature limit. The transition 
to this phase from the ordered state is obtained by pro- 
liferating the Zi vortices (which are line defects in the 
three dimensional classical model). In principle an ex- 
otic paramagnetic phase is also possible which would be 
the classical three dimensional analog of the spin liquid: 
this requires destroying the magnetic order without pro- 
liferating the Z 2 vortices. This paramagnet will be topo- 
logically ordered and will be separated from the trivial 
very high temperature paramagnet by a phase transition. 
In this classical context, the arguments of Ref. Q] apply 
to the transition between the Neel state and this topo- 
logically ordered paramagnet which will indeed be in the 
0(4) universality class. However neither this non-trivial 
paramagnet nor the corresponding transition were appar- 
ently accessed in the numerical calculations^. Equally 
it is hard to decide which of the two possible paramag- 
netic phases were accessed in the 2 + e calculations of Ref. 
I22I The distinction between the two phases is topological 
and hence sensitive to spatial dimension. This is hard to 
disentangle in the e expansion. 

Turning to the quantum problem at hand it differs 
from the classical stacked magnet in an important way. 
There are extra Berry phases that are sensitive to the mi- 
croscopic spin at each lattice site (i.e. spin- 1/2 or spin-1, 
etc). A close and familiar analogy is from the theory 
of collinear quantum antiferromagnets in two dimensions 
where such Berry phases spoil any general direct mapping 



to classical collinear magnets in one higher dimension. In 
that case in the semiclassical limit the Berry phases are 
associated entirely with singular topological configura- 
tions known as hedgehogs. A similar result also holds 
for the non-collinear quantum magnets of interest in this 
paper. The Berry phases are associated entirely with the 
topological Z 2 vortex configurations. As these Z2 vor- 
tices are gapped across the transition to the spin liquid 
the Berry phases play no role in the low energy universal 
critical physics. Indeed the 0(4) universality class will 
describe the Neel - spin liquid transition for all spin-S* 
magnets regardless of the value of S. On the other hand 
the phase obtained when the Neel state is disordered by 
proliferating the Z2 vortices will be strongly influenced 
by the Berry phases. For spin- 1/2 it is expected that 
the Berry phases will lead to broken lattice symmetries 
in the resulting paramagnel^. The nature of the tran- 
sition between such a VBS ordered paramagnet and the 
Neel state on the triangular lattice is not presently under- 
stood. Preliminary analysis suggests that an interesting 
"Landau-forbidden" second order transition may be pos- 
sible. However in the present paper we restrict ourselves 
to studying the transition to the spin liquid. 

IV. CONSEQUENCES OF CRITICAL THEORY 

A. Spin correlations 

We begin by considering the spin-spin correlation func- 
tion at the ordering wavevector which is readily accessed 
in neutron scattering experiments. This can be expressed 
in terms of the correlation functions of the vectors ni and 
n 2 

(S(r,r)-S(0,0)) = (m(r, r) • m(0, 0)) cos(Q ■ r). (15) 

We note that due to the symmetry of rotations between 
tii and n 2 , they will both have the same correlations. 
Close to the critical point, the correlation function has 
the usual power-law behavior 

(S(-q,o;).S(q,a;))~ ( ^_ fc 1 2)1 __ /2 , (16) 

where k = q — Q. Since ni and n 2 (and S) are composite 
operators in terms of the z a fields, one can expect that 
the spin-spin correlation function has a large anomalous 
exponent 77. Indeed, the large-TV calculation of 77 for an 
0(2N) theory givesi 



For the physical case N = 2, we have fj w 1.54. This 
is a large number compared to anomalous exponents of 
non-composite operators. 

To improve the large- N estimate of the anomalous ex- 
ponent, one can alternately perform Monte Carlo sim- 
ulations of the classical 0(4) nonlinear sigma model in 
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three dimensions since this model is in the same univer- 
sality class as the transition of interest. Ballesteros et 
alml measured the anomalous exponent t?t of the tenso- 
rial magnetization that is defined as 



M T = (VTiM 2 

Map = ^(0ia0i/3 - jSap), 



where 4>i a are four components of the 0(4) unity vector at 
the site i in three-dimensions. They found rj^ = 1.375(5). 
This result is directly related to our case because of the 
following argument. Using Eq. [5] and writing z-\ = (0i + 
i(f>2),zi — (03 + i(f>i), where Q are now real fields, we 
have 



ni = 



112 = 



-01- 

H02 
H03 



i + 0|- 
2(0104 



0203) 



(18) 
(19) 



We can see that the tensorial magnetization is related 
to ni and ri2 by symmetry and therefore the anomalous 
exponent of the tensorial magnetization is equal to the 
anomalous exponent of the vectors ni or 112. Then, from 
Eq. UH1 it is the same as the anomalous exponent of the 
correlation function at wavevector Q. 

We have confirmed this conclusion by computing the 
anomalous exponent of ni directly. We simulate the 0(4) 
nonlinear sigma model in three dimensions on the simple 
cubic lattice using a cluster algorithm—. We measure the 
anomalous exponent 77 in the following way. We calculate 
the total z component n\ — 2(0i03 — 0204) of the vector 
ni 



(20) 



fh and p — 7 = v(2—ff) for x), Fq is a universal function, 
t = \T — T c |, and a; is a universal exponent related to the 
leading irrelevant operator. Exactly at the critical point 
t = 0, one can write 



O = L p / u {a + cL- u ). 



(22) 



We assume that the correction to scaling is negligible for 
large lattice sizes and fit to a simplified scaling form 



O = aL p ' v . 



(23) 



We have performed simulations for large enough lattices 
(up to L max = 96). Fitting from L min = 20 to L max , we 
obtain fj = 1.373(3). 

The large value of fj is a reflection of the emergence of 
spinons as useful degrees of freedom at the critical point. 



B. Other competing orders 

The extra 0(4) symmetry actually has more striking 
consequences for the critical properties. It implies that 
operators other than the natural magnetic order param- 
eter will have enhanced power law correlators. Consider 
the magnetic order parameters n^a. As discussed in Sec- 
tion II I II they may be regarded as the first and second 
columns of a rotation matrix R. The third column of 
the rotation matrix is 113 = ni x 112. The enlarged 0(4) 
symmetry at the critical point implies that both left and 
right multiplications by an orthogonal matrix of this ro- 
tation matrix R are symmetries of the critical theory. 
Left multiplication is just physical spin rotation. How- 
ever symmetry under right multiplication implies that 
Hi, 2, 3 will all have the same correlations. In particular 
113 will have the same slow power law decay as 111,2 cal- 
culated in the previous subsection. We can express the 
vector 113 in terms of the spin operator as follows. Con- 
sider S r+ j x S r , where 1 is a unit lattice vector. Using 
and expanding the first spin operator around r, we have 



where the sum runs over all lattice sites. In our simula- 
tions, we measure the "magnetization" and "susceptibil- 
ity" 

rh = -^(Af). 

where L is the system size. We perform our simulations 
only at the critical temperature found in Ref. with 
very high accuracy, C = 1/T C = 0.935861(8). To find 
the critical exponents, we use finite size scaling analysis. 
Close to the critical point, observables satisfy the general 
scaling form 



O = LP' v Fa(tL x ' v ,L-"), 



(21) 



where O can be m or x, p is the scaling exponent corre- 
sponding to the operator O (p = — (3 = + fj)/2 for 



S r ^j x S r 



n 3 sin(Q • 1) 

derivative terms. (24) 



We can ignore the derivatives terms since their contri- 
bution to the correlation function is negligible. Consider 
the chirality vector 



CZ/j. — ^ ^ Sj. x Sj.' , 



(25) 



where the oriented sum is over an elementary plaquette. 
Using Eq. 1)24)1 . it is easy to show that 



C = 



sm ( 



2 sin — cos 
2 



n 3 cx n 3 . 



(26) 



Thus the vector spin-chirality correlation function 

(C(r,r) -C(0,0)) « (n 3 (r, r) • n 3 (0, 0)) 

= ( ni (r,r)-ni(0,0)). (27) 
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have the same power-law behavior at zero wavevector as 
those of spins at the wavevector Q. This highly non- 
trivial statement possibly provides a way to test the the- 
ory of Ref. Q in experiments and in model numerical cal- 
culations. 



C. Conserved quantities 

Some less dramatic effects of the enlarged 0(4) sym- 
metry are also worth noting. The 0(4) group has six 
generators and these will all be conserved at the critical 
fixed point. Three of these are just the conserved total 
spin Stot — J2i Si - they are the generators of left rota- 
tions of U. The remaining three are conserved only in 
the low energy critical theory (though not in the original 
microscopic model). They correspond to right rotations 
of U. Wc will denote these K a , a = 1,2,3. The K a 
transform as vectors under the group SU(2)r of right 
rotations of U . They generate rotations of n a amongst 
one another so that 



^^abc^c 



(28) 



As is well-known such conserved quantities have scaling 
dimension d S caie = 2 with no anomalous dimension 27 . 
Thus at the critical point at zero temperature 



(K a (r,r)K b (O,0)) 



S a b 



(| r |2 +T 2)2- 



(29) 



where r is the spatial coordinate and r is the (imag- 
inary) time coordinate. A similar result also holds 
for the conserved total spin. Away from the critical 
point at finite temperature in the quantum critical re- 
gion, the K a will continue to be approximately con- 
served. The exact conservation will be spoiled by irrel- 
evant operators that break the 0(4) symmetry down to 
SU(2) x lattice space group. Thus at T > in the quan- 
tum critical region we expect that the K a will diffuse 
upto a long length/time scale which will be determined 
by T and the distance to the critical point. It is possible 
that the presence of such extra nearly diffusive modes 
can also be looked for in experiments. 

But what do the K a correspond to in terms of the un- 
derlying spins? A useful guess for the answer is provided 
by examining the transformation properties of K a under 
the microscopic symmetries of the original lattice sys- 
tem. It is reasonable that any lattice operator that has 
the same transformation properties will have some over- 
lap in the continuum theory with the K a . With the fur- 
ther assumption that there are no other operators in the 
continuum theory with smaller scaling dimension that 
also have the same transformation, the K a will give the 
dominant correlations of such lattice operators. Similar 
considerations were also used in Ref. Il8l 

The transformation properties of K a under all physical 
symmetries is fixed by the commutation relation Eq. 1281 
The transformation properties of K a are related to the 



transformation properties of the matrix R. First, we list 
the latter ones. 



SU(2) 



spin 



R^OR, 



T T : R 



R 



cos(Q • 1) sin(Q • I) 
- sin(Q • 1) cos(Q • 1) 




R c ,1 : R — > R 



T : R^R 




where O is an 0(3) matrix and Tj is the lattice transla- 
tion along the lattice vector 1. We can use these trans- 
formations and the condition that K a are vectors under 
SU(2)k to obtain the transformations of K a 



577(2) 



where 



spin 

Rc-I 
T 



M 



(K U K 2 ,K 3 ) 
{K U K 2 ,K 3 ) 
(Ki,K 2 ,K 3 ) 
{K U K 2 ,K 3 ) 



{K X ,K 2 ,K 3 ), 
{K 1 ,K 2 ,K 3 )M, 
(K u -K 2 ,-K 3 ), 
{K X ,K 2 ,K 3 \ 



'cos(Q-l) -sin(Q-l) s 
sin(Q • 1) cos(Q • I) 
1 



(30) 



We can construct local spin operators that have the 
same symmetry properties as K a . We find that 

Ki(r) ~ cos(Q-r)S(r)-[S(r-x) + S(r + x)],(31) 
K 2 (r) - mn(Q-r)S(r).[S(r-*) + S(r + *)],(32) 
K 3 (r) ~ S(r - x) • [S(r) x S(r + *)]. (33) 

Thus K\ + iK 2 transforms identically with the Fourier 
transform of the bond energy at wavevector Q while K 3 
can be identified with the scalar spin-chirality along the 
chains. 



V. EXPERIMENTAL IMPLICATIONS 

A. Spin correlations 

Quite generally in the quantum critical region the spin 
correlations at the ordering wavevector must satisfy scal- 
ing - for instance as a function of lu/T. CS2C11CI4 orders 
at low temperature. Consequently data at the lowest 
temperatures and frequencies must be excluded from the 
scaling analysis. Furthermore it is precisely the very-low- 
T/low-w range that will also be most sensitive to the pres- 
ence of weak spin anisotropics (such as Dzyaloshinski- 
Moriya interactions). However the data at intermediate 
T and w will be less affected by such spin anisotropies or 
by the low energy long range order. 
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The value of fj found from Monte Carlo simulations 
can be compared to T\e estimated in experiments using 
a fitting procedure 5 . The experimental values fall in a 
range between 0.7 and 1. Thus the Monte-Carlo result is 
still larger than the experimentally estimated value. It is 
important to note the following caveat on this compari- 
son. The existing measurements of T\e were performed in 
inelastic neutron scattering experiments where both the 
wavenumber and frequency transfer to the sample were 
simultaneously varied. In particular typical data sets 
(such as for instance scan G of Ref. 5j) involve increas- 
ing the frequency transfer while the momentum transfer 
varies from points in the two-dimensional Brillouin zone 
far away from the ordering wavevector Q to points close 
to Q. In the context of the ideas explored in this paper 
the frequency dependence at fixed q far away from Q will 
not be dominated by the singular long distance critical 
fluctuations and will be highly non-universal. In con- 
trast the frequency dependence for q close to Q will be 
determined by the long distance critical fluctuations and 
will be universal. It is these latter fluctuations that are 
described by the calculated exponent fj. Thus the exist- 
ing measurements of f\E are quite possibly severely con- 
taminated by non-universal short distance effects. Hence 
the lack of quantitative agreement between fj and f}E is 
not surprising. However the qualitative observation of 
broad magnon linewidths is consistent with the large fj 
obtained in the theory. Future experiments will hopefully 
directly probe the fluctuations near the ordering wavevec- 
tor thereby allowing for quantitative comparison. 

As noted already in Ref. Q, NMR experiments may 
be a useful way to directly measure fj. For the nuclear 
relaxation rate we have 

7fr ~ T^- (34) 

Again this behavior will only obtain at intermediate tem- 
peratures. 

B. Detection of vector spin-chirality fluctuations 

A dramatic consequence of the critical theory studied 
in this paper is the enhanced vector spin-chirality corre- 
lations discussed in Section llV Bl How can they be de- 
tected in experiments? In this section we show how this 
may be done through polarized neutron scattering. Our 
proposal builds on and generalizes the pioneering ideas of 
Maleyev 28 in the context of classical non-collinear mag- 
nets. 

The part of the neutron scattering rate Rp that de- 
pends on the incoming neutron polarization Pj is given 
by 

J2p(k i ,P i -k / )~(P i .$)($.b), (35) 

where q = k/ — is the three- momentum transferred to 
the sample, and q is the corresponding unit vector. The 



vector b is given by 

b a = (36) 

where Sp~ f is the spin structure factor. Thus the polariza- 
tion dependent part probes the antisymmetric part of the 
structure factor. By the usual arguments it is clear that 
b a can be obtained from a calculation of the imaginary 
time Green function (here < r < (3 = 1/T) 

g(r,r) - (S(r,r) x S(0,0)), (37) 

g(q,»wn) = / e < * r -*"» T g(r,T). (38) 

J rr 

Specifically we have 



b(q,w) = 




FIG. 2: Dzyaloshinski-Moriya interaction is non-zero along 
oriented zigzag bonds that are denoted by vectors Si and ^2- 

In the context of this paper using Eq. [21 we find 

g(r,r) = (ni(r,r) x m(0,0))cos(Q-r) 
+ (n 2 (r,r) x m (0, 0)) sin(Q • r). 

Clearly such correlators are zero if there is full spin 
isotropy. In the specific case of CS2CUCI4 the presence of 
a weak Dzyaloshinski-Moriya(DM) interaction will lead 
to a non-zero result. It is known that the DM interac- 
tion is non-zero along oriented zigzag bonds. For a single 
triangular layer 

H D M = -D.53s P x(Sr+« 1 +S r+aa ). (40) 

r 

Here S1.2 are unit vectors along the oriented zigzag bonds 
as shown in Fig. [21 In addition Hdm is also staggered 
between different layers. The DM vector is oriented along 
the a axis and has magnitude D w 0.02meV « 0.05 J. In 
the continuum theory we may again use Eq. [21 to write 

Hdm « -d J d 2 r m x n 2 , (41) 

= d • / d 2 r n 3 . (42) 
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Here d oc D. Let us now try to evaluate g to linear order 
in d. We get 

g(x) = J tfV(ni(x) x ni (0) (d • n 3 (x'))) cos(Q • r) 

+ J d 3 x'((n 2 (x) x m(0)) (d ■ n 3 (x'))) sin(Q • r). 

Here x = (t,t),x' = (r',r') are space-time coordi- 
nates. The averages are evaluated in the isotropic the- 
ory. The first average vanishes due to the symmetry 
113 — > — 113, ni — > ni present in the action. We are there- 
fore left with 

g(s) = d a J d 3 x'(n 2 (x) x m{0)n 3a (x')) sin(Q • r). (43) 

Thus g is determined by the three point correlation func- 
tion of 1x1,2,3- Using spin rotation invariance we have 

g(r,r) = |sin(QT) J d 3 x'(n 2 (x)xn 1 (0)-n 3 {x')). (44) 

Such 3-point correlators of scaling fields are severely re- 
stricted by conformal invariance^. As n 12 ,3 all have the 
same scaling dimension A = (l+fj)/2, we have (at T = 0) 

(n 2 (x) x n x (0) • n 3 (x')) ~ ^7a^ _ X ,,A • ( 45 ) 

We now have to integrate over x' which is the coordi- 
nate of 113. It is easy to see that this integral is infra- 
red (IR) divergent. This is because at large x', the inte- 
grand behaves as -^k. With 2 A = 1 + 77 » 2.37, the 
integral over x' will diverge in the infra-red. Formally 
this just means that perturbation theory in the DM in- 
teraction diverges at T = and the correct answer will 
involve some fractional power of the DM vector. However 
for our purposes it is much more meaningful and simpler 
to go to finite temperature where this divergence will be 
cut-off. Before doing that it is useful to understand the 
origin of this divergence at T — better. 

The divergence comes from large x' 3> x. To discuss 
this limit let us keep x' fixed and bring x close to 0. 
Then we need the operator product expansion(OPE) of 
n2 x ni. This is a vector in spin space and the leading 
term will just be 113. Scaling requires the equation 

n 2 {x) x m(0) ~ ^n 3 (z). (46) 

There is only one power of A in the right hand side so 
that both sides will scale identically. Now if we calculate 
the correlator with 113(2;') we will reproduce the limit of 
the exact result Eq. above when x' is large. It is now 
clear that the infra-red divergence in the integral over 
x' actually is nothing but the divergence of the uniform 
susceptibility of the vector spin chirality at zero temper- 
ature. 

Armed with this insight let us discuss T > 0. Here 
we will assume that we are interested in the scattering 



at an external frequency u 3> T. This simplifies things 
because then the sole effect of T > is to cutoff the IR 
divergence without affecting the rest of the correlations. 

Using the OPE it is now clear that the answer is the 
Fourier transform of 

^Xcfc(T), (47) 

where Xch(T) is the uniform vector spin-chirality suscep- 
tibility at temperature T. From scaling we have 

Xoh(T) ~ (48) 

Using this we may straightforwardly calculate the inelas- 
tic scattering rate. For scattering right at the ordering 
wavevector Q and for lu T, we get 

R PM ~ CJ 3-Ar.3-2A (49) 

with 2A = 1 + fj w 2.37. 

This is a definite prediction that can possibly be tested. 
In the actual experiments it will be necessary to take into 
account the staggering of the DM interaction between the 
different layers. Thus it is best to choose the a component 
of the 3-vector q to be it. 



VI. SUMMARY AND CONCLUSIONS 

In this paper we have pursued a concrete version of the 
idea that CS2CUCI4, though magnetically ordered at low 
temperature, may nevertheless be proximate to a spin 
liquid phase. Such proximity suggests that the inter- 
mediate length/time scale physics of CS2CUCI4 may be 
governed by a quantum critical point between a magnetic 
spiral and a genuine two dimensional quantum spin liq- 
uid. A theory for such a quantum phase transition was 
obtained in Ref. [lj for a simple (Z 2 ) spin liquid with 
gapped bosonic spinons. We showed that the spin cor- 
relations at the quantum critical point are characterized 
by a large anomalous exponent fj 1.37. This is qualita- 
tively consistent with the broad power law tails observed 
in inelastic neutron scattering in CS2CUCI4 . Further we 
showed that the enlarged 0(4) symmetry^ at the critical 
fixed point has some remarkable consequences. The vec- 
tor spin-chirality has the same slow power law decay as 
the natural magnetic order parameter. This sharp qual- 
itative observation should be of great use in confirming 
(or ruling out) the applicability of the theory of Ref. 1 
to Cs 2 CuCl4 . Building on Ref. l28l we showed how 
polarized inelastic neutron scattering can be used to di- 
rectly detect the vector spin-chirality correlations. It is 
our hope that future experiments will be able to use these 
results to clarify the physics behind the interesting prop- 
erties of Cs 2 CuCl4 . 
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